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Abstract 

This is a review paper on recent work about the connections between rough path theory, the 
Connes-Kreimer Hopf algebra on rooted trees and the analysis of finite and infinite dimensional 
differential equation. We try to explain and motivate the theory of rough paths introduced by 
T. Lyons in the context of differential equations in presence of irregular noises. We show how it 
is used in an abstract algebraic approach to the definition of integrals over paths which involves 
a cochain complex of finite increments. In the context of such abstract integration theories we 
outline a connection with the combinatorics of rooted trees. As interesting examples where 
these ideas apply we present two infinite dimensional dynamical systems: the Navier-Stokes 
equation and the Korteweg-de-Vries equation. 



Keywords: Rough Path theory, Connes-Kreimer Hopf algebra, driven differential equa- 
tions, Navier-Stokes equation, Korteweg-de-Vries equation. 

1 Introduction 

Rough path theory has been developed by T. Lyons for the analysis of the map $ : a; i-^ y which 
sends a vector- valued driving signal x to the solution of the differential equation 



dyt ^ X! f^ivt)^^^, yo^y (i) 



where y is the inital condition. This equation has a well defined meaning as a non-autonomous 
ODE when xt is differentiable in t. In applications however it is interesting to consider more general 
driving signals, e.g. when taking into account random perturbations of dynamical systems where 
Xt is the Brownian motion which is almost surely non-differentiable. Indeed the application to 
stochastic differential equations (SDEs) has been the main motivation in the developments of the 
theory. The standard approach is then to understand eq. ([l]) as the integral equation 



yt = y- 



J2faiys)dx:, (2) 
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provide a well-defined meaning for the stochastic integral in the r.h.s. and the proceed to solve the 
equation by standard fixed-point methods. Usually the integral above can be an Ito integral or a 
Stratonovich one, but other more exotic choices are possible in other stochastic contexts (e.g. the 
Skorohod integral, the "normal ordered" integral, etc. . . ). Whenever we speak of a solution y to 
the SDE ll]) we mean suitable random fmiction y for which f{y) can be integrated against dx in 
an appropriate sense. 

Lyons' basic observation was that, in the case of a smooth control x, the solution of the differ- 
ential equation is a well-behaved function of the iterated integrals X of x: 

Xfs^ ■■■ dxZ---dxll (3) 

J s J s J s 

where we denote with a the multi-index (oi, . . . ,a„) with \a\ — n its length. In some sense the 
iterated integrals encode the local behaviour of the path x well enough to faithfully recover its 
effect on the solution y. In a system perspective, the iterated itegrals provide a canonical set of 
coordinates for the analysis of non-linear systems much like the standard Fourier coefficients are 
natural coordinates for linear ones. 

Like Fourier coefficients, iterated integrals enjoy nice relations upon concatenation of paths: 
given three times s < u < t we have the celebrated Chen relations ^Ij between iterated integrals: 

n-l 

xf--" = x^r"'' + Kr''- + E ^°u"-"'=^:r"-""- (4) 

fe=i 

These non-linear equations plays a fundamental role in the development of the theory. 

The plan of the paper is the following: in Sect. [2] we motivate the simplest instance of a rough 
path in the context of integration against smooth approximations of an irregular path. Sect. [3] 
we describe an algebraic approach to integration and its use to define integrals against a rough 
path. This approach is not standard from the point of view of rough paths. More conventional 
expositions of the theory can be found in [21 [3] and [4]. In Sect. [4] we describe the combinatorics 
of iterated integrals of very general type via the Connes-Kreimer Hopf algebra on rooted trees 
following ^ . This will ultimately allow to define and solve differential equation associated to such 
integrals. Finally in Sect.[5]we apply the objects and the related combinatorics of algebraic integrals 
to the study of infinite dimensional differential equations via two prototypical examples: the 3d 
Navier-Stokes equation as studied in [6l and the Id periodic Korteweg-de-Vries equation [7j. 

2 Renormalizable theories of integration 

The analysis of the integral equation ([2]) in the case when a; is a non-differentiable function can be 
split in two parts: the definition of the integral on the r.h.s. and the fix-point argument. Let us 
concentrate on the first part. Assume given a continuous function x : [0, T] — R" which we assume 
only 7-II61der continuous, i.e. for which the following estimate holds 

\xt-x,\<C\t- sW i,se[0,T] 

(where the least constant C define the Holder norm |ja;||-y) and for simplicity let us restrict in this 
section to the case 1/3 < 7 < 1. Let us pose the problem of giving a "reasonable" definition of the 
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integral 



for some smooth one- form (/? : R" — > R". Proceeding by approximation we consider a family 
{a;(e)}£>o of smooth approximations to the path x and let 



IVp\{^)ts = / ^p{x{e)u)dx{£)u- 

J s 



It is easy to convince ourselves that in general we do not have any control of these integrals if 
x{e) X when e ^ in the "/-Holder norm. It is then a remarkable fact that all the possible 
integrals obtained varying the function all converge at the same time (provided is sufficiently 
smooth) when the approximated iterated integral of order two: 



J s J s 



converges to a function X : [0, T]^ M" ® R" in the sense that 

sup i— — I ^0 as £ ^ 0. 

0<s<t<T |t-s|^T 

So in some sense we can claim that the integration theory agains the path x is well defined as long 
as we are able to control the convergence of the iterated integral of order two X°''^°'^: all the other 
integrals (and more fundamentally also all the higher order iterated integrals) will turn out to be 
nice functions of the data given by the path x and X"-^'''^ . 

It is suggestive to understand this phenomenon as a very simple example of "renormalizable 
theory" where all the quantities of interests have well defined meanings as functionals of a finite 
number of fundamental objects whose intrinsic determination remains outside the scope of the 
theory itself. 

In this limited context we face the appeareance of the simplest non-trivial example of a rough 
path: the couple (X", X°-^°-^) where X^^ = a;^ — is a 7-rough path, i.e. a path and some additional 
information in the form of "iterated integrals", for which a complete theory of integration and 
differential equations can be constructed (as we will see shortly). The fact that it is enough to 
consider only the second order integral is due to our hypothesis that 7 > 1/3. 

Note that once the limit has been taken the object X"-^"-^ is no more an iterated integral (in a 
classical sense) and can be characterized more abstractly by the following two properties 

1. the Chen relation: 

sup p^<oo. (6) 

and more interestingly there could be more than one possible choice for this object compatible with 
these two conditions leading to different integration theories. 



2. a regularity condition 
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3 Algebraic integration 



To understand how the iterated integrals X conies into play in the definition of I[(p\ we need some 
tool which allows us to analyze the "local" (with respect to the parameter) behaviour of the integral. 
Specifically we want to expand the integral in a short interval as 

I[ip\ts = '^>a{.Xs)Xl + dbipa{Xs)Xtt + Us (7) 

where r stands for some remainder term which we hope will be of higher order in |i — s| than the 
other terms. The main property of the integral I[(f\ is its trivial behaviour under splitting of the 
integration interval: I['^]ts = I[^]tu-'rl[^]us for s < u < t, this of course means also that the integral 
can be written as the increment of a function ts = ft — fs (for example taking ft — /[(p]ot). In [8J 
we introduced a cochain complex {C■^,,5) which encodes this basic property of integrals. For n > 1 
let Cn be set of functions a G C([0, r]";R) such that at^---t„ = if = i^+i for some 1 < i < n — 1. 
Elements in C„ will be called n-increments. Define a coboundary d : Cn Cn+i as 6ats = at — as 
for a S Ci, Satus — o-ts — o-tu — o-us for a E C2 and so on. Let ZCn = KerS fl C„ and BCn = ImS fl C„. 
It is easy then to verify that 66 = and that the complex {C^,6) is exact, i.e. that lm6 — Ker6 at 
any C„, n > 1. In particular a 2-increment a G C2 is the increment of a function / G Ci if and only 
if 6a = 0. So at C2 the coboundary measures the degree of "exactness" of 2-increments. Moreover 
a key fact is the following: if 6a is suitably small then there exists only one function / (modulo 
constants) such that 

6f^a + r (8) 

where the remainder is small. To be more precise we need to introduce the relevant notion of 
"smallness" . We say that a G if 

.. \ats\ 

sup ■; ;— < 00. 



And similarly 6 G C3 if 



III \btus\ 

\bh ■■= sup < C30. 

s<u<t \t — S\' 



Both C2 and are Banach spaces when endowed with the norms || • Define = Uz>iC^. More- 
over if g G Cn and h G Cm then we write gh for the element of Cn+m-i such that {gh)tj^.--t„^^_i — 
9ti---t„ht„---t„+^_i- 

The key result in this theory of increments is the existence of a sewing map which provide a 
natural inverse operation to 5: 

1+ 



Theorem 3.1 (The sewing map). There exists a unique bounded and linear map A : ZC^ — > 
such that 6 Ah = h for any h G ZC^ 



1+ 
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In particular the map (1 — A6) : C2 — > C2 projects (in essentially a unique way) suitable elements 
of C2 to ZC2 — BCi. the space of increments of functions. Then going back to the decomposition ([8]) 
and assuming that 6a G Cg^ we can form the exact 2-increment a — K6a = 6f and obtain that the 
remainder r is given by K6a G C2^ . It is also easy to see that, the decomposition ([5]) is unique if we 
require r G Cj^, indeed if two such decomposition exists, (/, r) and (/',r'), their difference satisfy 
6{f — /') = r — r' G C| for some z > 1 and f — f would be a function of z-H61der class ans since 
z > 1 that this function would take the constant value zero. 
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With this notions at hand we reahze that eq. ([7]) is nothing more than an instance of a decom- 
position similar to ([8]). As we have aheady seen we can determine both /[y;] and r at once using 
only the well-defined remaining terms in the r.h.s., so letting ats = (paixs)^^^ + db(pa{xs)Xf^ , in 
order to be able to apply A, we need to require that Sa G Cl~^ . By some easy computation using 
the fact that 6 satisfy some kind of Leibniz rule, we get 

Satus = ~Sipaix)tuX^^ - Sdbipaix)tuX^^ + dbipaiXs)SX^^^ 

If we now exploit the Chen relation for A"*"" we can simplify this expression further 

Satus = -[SiPaix)tu - dbipa{Xs)SXl]X^^ - Sdbipa{x)tuX^: (9) 

Assuming that (p is sufficiently smooth (C^ is enough) we obtain easily that 

Scpaix) - dbipa{^)SX' € Cl\ Sdbifiaix) G C] 

SO that taking into account all the regularities we end up with Sa € C^'^ . Then our assumption on 
7 ensure that 87 > 1 and we can prove that there exists a unique couple (/, r) with f £ Ci and 
r e such that 

6f = ipaiXs)X^^ + dbipai^^)Xt + rts- (10) 

By construction / depends only on the 7-rough path X and on ip. We can then define the integral 

by 

I[ip] =Sf={l- AJ)[XVa(x) + X^^db^x)]. 

To motivate the fact that the decomposition (|10p is a sort of renormalized integral we can make 
the following observation. Take a partition {ti} of [s, t] of size A and consider the Riemman sums 
Sa = X^i v(2;ti)^ti+iti- In general we have no mean to say that these sums are convergent as 
A — > 0. However given a 7-rough path X we can perform a subtraction to these sums and define 

S'a^T. [p-i='u)K^,u+dbV-ixu)X^^^.u 

i 

then using the decQniDOsitioi JTUl as A ^ we have the hmit 

i i i 

since the first sum telescopes and the second is easily show to converge to zero. 

Moreover the regular dependence on the data ensure that smooth approximations I[ip]{e) con- 
verge to I[ip] as here defined as long as we can prove the convergence of approximating path (and 
its second order iterated integrals) to the rough path X. 

Exploiting the sewing map, the integral equation ([2]) in presence of a 7-rough path with 7 > 1/3 
can be understood as a fixed-point equation for an unknown y E Ci: 

Sy^ii- AS)[x-My) + x^^dbMy)^%)] 

these equation can be solved by a standard iteration method in a suitable subspace of Ci ([5] for 
details) . 

It is also possible to construct higher order iterated integrals starting from low order ones where 
for iterated integrals we means simply object which obey Chen's relations Take for example 
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the third order object X^^"^"^ G C2, eq. ^ can be written as a statement about the coboundary 

of X'^^"'^'^^: 

and in the above hypothesis on X it is easy to check that the r.h.s. belongs to ZC2"' C ZC2^ so 
that it is in the domain of A and we can define 

as the unique solution in C^^ of the Chen relation. Iteratively this allows to construct all the higher 
order iterated integrals. Some more work allows to prove a uniform estimate on the growth of the 
norms involved in the procedure [21 [8] : 

m\rah<Ci^ (11) 

for any multiindex a. 



4 Trees 

What happens if 7 < 1 /3? The obstuction to the exactness of the increment a in eq. (O will no 
more belong to the domain of the sewing map. Indeed in eq. (O we cannot anymore expect that 
the remainder belongs to and at the very least this would affect our argument for uniqueness. 
We are then forced to proceed and expand further the integral, or from another point of view, to 
add some counterterms to remove large contributions to r. 

To understand the general structures that we need if we proceed further in the expansions it is 
better to take as working bench the more difficult case of the differential equation The series 
solution will be indexed by rooted trees: a phenomenon which is present already for solutions to 
the ODE dy/dt = f{y) [i[10]. 



4.1 Labelled rooted trees 

Given a finite set C, the /^-labeled rooted trees are the finite graphs where labels of C are attached 
to each vertex and where there is a special vertex called root such that there is a unique path from 
the root to any other vertex. 

Some examples of rooted trees labeled by £ = {1, 2, 3} are 

« K V.' if' 

Trees does not distinguish the order of branches to each vertex. Given k £-decorated rooted trees 
Ti, • • • , Tfc and a label a E C we define t — [ti, ■ ■ ■ , Tk]a as the tree obtained by attaching the k 
roots of Ti, • ■ • , Tfc to a new vertex with label a which will be the root of r. Any decorated rooted 
tree can be constructed using the simple decorated tree •a (a G ^) and the operation [■••], e.g. 
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Denote Tc the set of all L decorated rooted trees and let T the set of rooted trees without 
decoration (i.e. for which the set of labels L is made of a single element). There is a canonical 
map 7c — > T which simply forget all the labels and every function on T can be extended, using 
this map to a function on Tc for any set of labels L. Let | • | : T — > M the weight which counts the 
number of vertices of the (undecorated) tree defined recursively as 

|.| = 1, |[Tl,...,Tfc]| = l + |Tl| + ...+ |rfe| 

moreover the factorial 7(t) of a tree r is 

7(») = 1, 7([ti, . . . ,Tfc]) = |[ri, . . . ,rfe]|7(Ti) • ■ •7(rfe) 
finally the symmetry factor a is given by the recursive formula (t(t) = 1 for |r| = 1 and 

S[t\ . . . ,t'=) 

where S{t^ , ■ ■ ■ , t*^) counts the number of different ordered /c-uples (r^, • • • , r'') which corresponds 
to the same (unordered) collection {r^, • • • ,t'^} of subtrees. The factor kl/S{T^, . . . ,t'°) counts the 
order of the subgroup of permutations of k elements which does not change the ordered fc-uple 
(r^, • • • , r*^). Then (t(t) is is the order of the subgroup of permutations on the vertex of the tree 
T which do not change the tree (taking into account also the labels). Another equivalent recursive 
definition for a is 

fT([(ri)"i ■ • • irT'U) - ni! • • • nMrT ' ' ' ^(^fe)"'= 

where , . . . are distinct subtrees and ni, . . . , the respective multiplicities. 

The algebra ATc is the commutative free polynomial algebra generated by {1} U Tc over M, 
i.e. elements of ATc are finite linear combination with coefficients in M of formal monomials in the 
form T1T2 ■ ■ ■ Tn with ri, . . . , t„ G or of the miit 1 G ATc- The tree monomials are called forests 
and are collectively denoted Tc ; we include the empty forest 1 g Tc ■ The algebra ATc is endowed 
with a graduation g given by giji ■ ■ • Tn) = |ti| + • • • + |r„| and g{l) = 0. This graduation induces 
a corresponding filtration of ATc in finite dimensional linear subspaces AnTc generated by the set 
^2 of forests of degree < n. 

Any map f : Tc ^ A where A is some commutative algebra, can be extended in a unique way 
to a homomorphism / : ATc ^ A by setting: /(n • • ■ t„) = /(ri)/(T2) • • ■ /(r„). 

In the following we will use letters r, p, ct, . . . to denote trees in Tc or forests in J-c, the degree 
g{T) of a forest r G J-c will also be written as |t|. Roman letters a,b, c, ■ ■ ■ ^ C will denote vector 
indexes (i.e. labels) while a,b, . . . will denote multi-indexes with values in C: a = (ai, . . . , a„) G 
with \a\ = n the size of this multi-index. 

4.2 Iterated integrals and the Connes-Kreimer Hopf algebra 

Given a smooth path x G C^([0, r],R") we can canonically associate to it a family of 2-increments 
A"^ indexed by trees labelled hy C — {1, . . . , n} hy the iterated integrals 

X*; =xt- X,, = / X2--- Xl-dxl (13) 

J s 

Iterated integrals like ^ corresponds to "hnear" trees t = [[• • • [•a„]a„_i • ■ •]a2]ai and by abuse 
of notation we will continue to write X"! '"" for such X'^ . The generalization of the Chen's 
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relations ^ involves the Hopf algebra structure on on labelled trees essentially introduced by 
Connes and Kreimer which we now brefly describe. 

On the algebra ATc we can define a counit e : ATc — > R as an algebra homomorphism such 
that ~ 1 and e(T) — otherwise and a coproduct A : ATc ATc ® ATc as the algebra 
homomorphism such that 

A(t) -l®T + ^(B;®id)[A(B^(r))] (14) 

on trees t &Tc, where = »a and S+(ti • • • t„) = [ti - ■ ■Tn]a and Bt is the inverse of or 

is equal to zero if the tree root does not have label a, i.e. 



B<L[BX{t,---t^)) 



Ti ■ • • T„ if a = & 
otherwise. 



Endowed with e and A the algebra ATc become a bialgebra, there exists also an antipode S which 
complete the definition of the Hopf algebra structure on Tc as described by Connes-Kreimer [11] 
(in the unlabeled case). 

We will often use Sweedler's notation for the coproduct At = ^ T(i) ®T(2) but we also introduce 
a counting function c : Tc y< Tc y. Tc N such that At = J^peTc a<^Tc ® ^ moreover 

it will be useful to consider also the reduced coproduct A't — At — 1®t — t®1 with counting 
function c'. 

The generalization of eq. ^ reads as follows: 
Theorem 4.1 (Tree multiplicative property). The map X satisfy the algebraic relation 

SX-" = X^'^-^l creATc (15) 

Let us give an example in one dimension (c? — 1) where trees are not decorated. The forests 
with |t| < 3 are 

The action of the reduced coproduct on these forests and the corresponding action of the coboundary 
on the iterated integrals are given by 



A'l = . (g) 
A'(..) = 2. ( 

A'J = I ®.+.«>I 5^us = ^lu^us + ^tu^us 

A'(.i) sxll = xi^xi:+xt:xi,+xlxi,+xi^xl 

A'(.^)=3.^(».+3.«i.=' _ ov*^ V* I ov* Y-*^ 

A'y =»«i»»+2i(»» V t 

Ay* — -U 9 

As a further example consider the iterated integrals T'^ associated to the identity path 

[ ] /"* 
T^g — t s, T^g — I T^tl ■ • ■ T^,^du 



us 
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By induction it is not difficult to prove that T^^ = (t ^ s)I'^I(t!) ^, so applying Thm. 14.11 to T'^ we 
get a remarkable binomial-likc formula for the Conncs-Krcimcr coproduct 

Iterated integrals of this sort appears naturally when trying to expand in series the solution 
of driven differential equations: for any analytic vectorfield / : M'' ^ M" and any smooth path 
X G C"'^([0, r],R"), the solution of the differential equation dyt — 'Ylaec fo-^V^^'^^t i with = r/ 
admit locally the series representation 

Syts = E -^<l^^ir){ys)x;^, yo = v (17) 

where we recursively define functions (f)^ {t) as 

fc 

0^(-a)(e) -/a(0, ^f{[r'---rXm - faM...H{i)X{W{r'mP. 

beICi:\b\ = k «=1 

Note that using the regularity of the path x the iterated integrals X"^ can be simplified: indeed 
Chen [IJ proved that products of iterated integrals can be always expressed as linear combination 
of iterated integrals via the shuffle product: 

xa.-a^xir'-^ J2 x:r'-^- (is) 

c<£Sh(a,b) 

where given two multi-indexes a = (ai, • • • , a„) and b — (bi, ■ ■ ■ , bn) their shuffles Sh(a, b) is the set 
of all the possible permutations of the (n + TO)-uple (ai, . . . , a„, 6i, . . . , bm) which does not change 
the ordering of the two subsets a, b. Using relation (fTSj) every X'^ can be reduced to a linear 
combination of standard iterated integrals {X°}a:. 

It is however interesting that many costructions related to integrals and solution of driven 
differential equations do not depend on eq. (|18p being valid. For example in the theory of the Ito 
integral eq. (|18p does not hold and more complex algebraic relations have to be considered: if x is 
a multidimensional Brownian motion and X defined via Ito integration we have 

X^Xl = ^tt + ^t'r + 5ab{t - s). 

Let us clarify the algebraic framework in which we consider possible integration theories. It 
turns out that the only data we really need to build a family {X'^}r£Tc satisfying ([T5|) is given by 
a set of linear maps {/"}ag£ defined on a domain Dj C (the unital algebra obtained extending 
C2 with a unit e and considering the point-wise product) to C2 satisfying two properties: 

1. I{hf)ts = Iih)tsfs for aU h eT>i,f <eCi and where {hf)ts = hsfs] 
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2. 61{h)tus = I{e)tuKs + Iih^'%uh:tl when h^Vi with 5htus - E» ^tn and h^ ' e P/ . 

Then given a family {/"}ae£ of such integral maps on a common algebra I? C C2 we can associate 
to them a family {X'^^T-eJ^c. recursively as 

r _ 1 _ fc T It- ^ h 1 t 

^ts lejts, ^ts — 1 {J^ }ts, ^ts — ^ts ■ ■ ■ ^ts ■ 

In this way we estabilish an algebra homomorphism from ATc to a subalgebra of C2 generated by 
the X'^-s. This homomorphism send the operation B°j_ on ATc to the integral map I'^ on C2 in such 
a way that Theorem 14. II holds. 

Given X'^ for any t ^ Tc with |t| < n and such that < C for some constant C and 

7 e (!/("- + 1), 1) we are able to extend to any r with |r| > n the "integrals" X'^ e C2'^' in a 
unique manner solving the equation SX'^ = X^ with the aid of the sewing map. This extension 
will satisfy the bound 

II^^IUi.i < c,d;^q,{T) (19) 

for any r with Ci, C2 two finite constants and q-^ir) a function satisfying 

The actual asympotic behavior of for |t| ^ 00 is not yet know but we conjecture that it should 
hold the equivalence 

q^ir) ^ (r!)-^ (20) 

(see [S]). This asymptotic behavior is satisfied on the subset of linear trees and is consistent with 
our results on tree-indexed iterated integrals in the context of 3d Navier-Stokes equation ^6^ (see 
also Sect. [5]). 

We call the object {X'^}r satisfying and pop for some 7 £ (0, 1) a {-j-) branched rough path. 
By the above consideration it is clear that it is determined by the finite subset {X'^}^|t-|<i. 

Assuming that the integrals {/°}aG£ generate a 7-branched rough path in the above sense, then 
we can extend the integral to a larger class of integrands. Take a set of constants {/io}r and consider 
the path h G Ci defined as ht — J2t ^o^Mj then 

shts = E E ^tl^'^ir = E ^'(^' '^Woxp.x!, = E ^^s^ts (21) 

where we introduced paths £ Ci by hi — ^ ^ c'(p, r, a)h^Xgf^. We must have 

= d^X^h'' = E^^*''^'"'"^^ - X''6h^ 

T T 

SO we can check that 

6h'' = E c'(/0, T, ct)X''/i'' (22) 

p,r,cj 

also holds. At this level these relations are formal since they require an infinite number of terms 
to hold, however exploiting our analysis of the (C*,(5) complex and the sewing map we can work 
modulo Cj"*" and neglect the terms in the expansions which involve X'^-s with 7|t| > 1. A more 
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accurate analysis (which can be found in ^) reveals that for the purpose of integrating and solving 
differential equations we can actually work modulo larger terms but for the sake of clarity we refrain 
to do this here. 

We call a controlled path (by X) any path h which satisfy modulo the eq. ([2T|) together 
with eq. ([22]) for all trees r with 7|t| < 1. The integral / can then be extended to integrate any 
controlled path. Indeed since he = eh + Sh it is a consistent definition to set 

/"(/i) = r(/ie) = r(e)/i + r(5/i) = ^ r(x^)/i^ + r(7^) 

7|r|<l 

7|t|<1 

where we exploited the properties of the integral /° and the definition of X and where TZ G C^^ 
is a generic remainder term (possibly different from line to line). Controlled paths are also stable 
under mapping by sufficiently smooth functions / (the degree of differentiability depends on 7): i.e. 
if y is a controlled path, then z = f{y) is again a controlled path with an explicit formula for the 
coefficients z'^ in term of the derivatives of / and of the coefficients y'^ . Using these properties we 
can consider the differential equation for the integrals / and vectorfield /: 

in the space of controlled paths y and solve it by a fixed-point argument. 

5 Infinite dimensional dynamical systems 

In this section we would like to show how the ideas and the tools described before could be applied 
in the context of infinite dimensional dynamical system by introducing operator-valued iterated 
integrals (and rough paths) which by their non-commutative nature are intrinsically indexed by 
trees. In particular we will discuss two different examples: the Id periodic Korteweg-de-Vries 
(KdV) equation and the 3d Navier Stokes (NS) equation. In the first case we will exploit the 
increment complex to analyze perturbatively the solution for irregular initial data , in the second 
case we will use series expansion over trees of the solution to analyze the dynamics for large times . 

As before we will restrict ourself to an overview of the results and to a sketch of the arguments. 
The interested readers can find a rigorous discussion elsewhere [71 |B] and he can refer to [H] for 
extension of these consideration covering the analysis of stochastic partial differential equations 
driven by irregular noises. 

5.1 The KdV equation 

The Id periodic KdV equation is the partial differential equation 

dtu{t,0+dlu{t,0 + \d^uit,0^ ^0, w(0,O = "o(e), (t,e)eKxT (23) 

on the torus T — [— tt, tt]. When the initial data is not smooth this equation must be interpreted in 
the integral form 

1 /■* 

u{t) = Uit)uo + - U{t- s)d[u{s)]^ds 

2 Jo 
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where U is the "free propagator" given by the hnear equation dtU{t) = d^U{t), d denoting the 
spatial derivative, and where we consider the solution m as a path in a space of functions. By going 
to the interaction picture vt — U {~t)u{t) we get 

1 f* 

Vt = uo + - U{~s)d[U{s)d,]^ds 
^ Jo 

Then the Fourier coefficients {vt{k)}kez of vt satisfy the equation 

Vt{k)^Voik)+'-^ f e-'^''"-''"-'"^'vs{k,)vs{k2)ds, te[0,T],k^O. (24) 

k=ki+k2,k„^0''^ 

We restrict our attention to initial conditions such that wo(0) = 0. By calling X the bilinear 
operator in the r.h.s. this equation takes the abstract form 

ft 



Vt = Vs + / Xa{Va,Vcr)da, t, S e [0,T]. 

J s 



where the paths take values in the Hilbert space Ha of complex Fourier coefficients (p{k) with 
(^(0) = and (p{—k) = (p{k) endowed with the scalar product {(pi,<f2)a = Sfc^o ^^"''^l(~^)'''2(fc)■ 
By iteratively substituting the unknown in this integral equation we obtain an expansion whose 
first terms looks like 



Vt ^ Vs + J daX^{vs,Vs) + 2 J daX^rivs, J daiXai{vs,Vs)) 

nt per n(T 

+ / daXai d(JiX„^{vs,Vs), da2X^^{vs,Vs)) (25) 



+ 4y daXaivs, J daiXaiivs, J da2X^^{vs,Vs)) + 

where rts stands for the remaining terms in the expansion. Denote with BT the set of (unlabeled) 
planar rooted trees with at most two branches at each node. A planar tree is a rooted tree endowed 
with an ordering of the branches at each node. Then each of the terms in this expansion can be 
associated to a tree in BT and we can define recursively multi-linear operators X'^ as 

X^^{ifl,ip2) = X„{Lpi,Lp2)d(j; 

and ^ 

Ml ^ '('l^'l' • ■ ■ ,<Pm+n) = XcriX^g{(pi, . . . , (fm) , i'fim+l , ■ ■ ■,ipm+n))da. 

Eq. (Uni) has the form 

5vts = X'iv^'^s + Xliv^'^s + xliv^'^s + X^iv^^^s + rts (26) 



12 



as an equation in C2 where increments take values in Ha and where we let w^" — {vs, ■ ■ ■ ,Vs) {n 
times). The operators X'^ satisfy multiplicative relations 

dXi{ipi,Lp2,ip3) = X'{X'{ipi,tp2),ip3), 
SX^{lfl,ip2,f3,f'l) = X'{X^{ipi,ip2,'^3),'P4) + X^{X'i(pi,ip2),'P3,'Pi), 

and 

6X^{ipi,ip2,(P3,(P4) = X'{X'{ipi,Lp2),X*{ip3,(f4)) 

+ X^i'fl,V>2,X'{(f3, (^4)) + Xi{(p3, (^4, X'{lPi, LP2)) 

where we used the symmetry of the operator X to obtain this last equation. These relations have 
much in common with the analogous relations for branched rough paths, however here the additional 
information of the position of the various arguments must be taken into account in the combinatorics 
of the reduced coproduct. It would be interesting to determine a Hopf algebra structure on BT 
which could account this additional information in a general way. 

Note that X is unbounded on Ha while it is possible to prove that X' and xl are bounded, in 
particular for the simplest of them, X* , we have the bound 

\X's{'Pl,f2)\a < C\t - Sp\(pi\a\(p2\a 

with the parameters 7 and a satisfying 

7 < 1/2, a e M with 7 - a< 1 and 87 - a < 3/2. (27) 
Moreover we are able to decompose the second member of the hierarchy xl as 

Xi{(pi,ip2, ^3) = X^{^l,V2,^3) + '^{f 1,^2)^3 + '^{^l,Vi)^2 (28) 

where for the same range of 7, a as in (|27|) we have 

\x}^{^l,tp2,^3)\a < C\t - S\^^\ipi\a\^2\a\^3\a 

while the operator ^ts ■ Ha x Ha — > C is bounded only for a > —1/2 and 

\^ts{Vl,iP2)\ < C\t - s||(/7l|a|'^2|Q- 

The analysis of the higher order operators has not yet been performed. However already at this 
stage something can be said if we take 87 > 1 since we are naturally led to consider eq. (pS)) as a 
increment equation and rewrite it using the sewing map and the operators X (only up to second 
order) obtaining the equation 

5w = (l-A(5)[X*(t;^2)^ ^1(^x3)] (29) 

which can be solved by fixed point methods in Ha for any a > —1/2 (cfr. (|77|) and the condition 
on $). 

The condition a > —1/2 is essentially imposed by the operator $ appearing in the decomposi- 
tion (HSl) of xl. This constraint it is linked with a resonance phenomenon which apperars in the 
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scattering by the non-Unear term involving four waves and which hint to the fact that the KdV 
equation is not uniformly wellposed in Ha for a < —1/2 [l3[fT4] . 

We have replaced the differential and integral approach to the study of this equation by an 
approach based on an operator valued rough path and the increment complex. It is interesting 
then to look how the properties of the dynamical system reflect in this unusual approach. As an 
example let us consider conservation laws. 

The KdV equation formally conserves the Hq norm. This conservation law imposes additional 
algebraic relations to the operators X: it is not difficult to prove that we have 

{fl,Xts{V2,f3))o + {ip2,X^g{ifi,ip3))o + {ip3,X^g{ip2,ipi))o = (30) 

and that 

2{ip,XUip,ip,^))o + {Xtsiip,ip),XUip,ip))o=0, (31) 

where all the test functions belong to Hq. To see that these two relations imply the Hq conservation 
law for solutions we will prove that 5{v,v)o = when v satisfy 

Sv = X'{v''^)+xl{v''^)+n (32) 

cfr. eq. dlH). Let us compute explicitly [6{v,v)o]ts = {vt,vt)o - {vs,Vs)q = 2{Svts,Vs)o + {Svts,Svts)o 
. Substituting in this expression the equation (j32p we get 

[S{v, v)o]ts = 2{X'g{Vs,Vs) + x}g{Vs,Vs,Vs), Vs)o 

+ (Xlg{v,,v,),X^g{v,,v,))o+n. 

The relation ([50]) implies that {vg, X'g{vs,Vs))o = while eq. (pij) allows to cancel the X^ term 
with the quadratic X* term. After the cancellations the increment of the Hq norm squared is then 
[(5(w,v)o] S but this means that it must be zero and that {vt,Vt)o — {vo,vo)o for any i > 0. 

5.2 Navier-Stokes-like equations 

We consider the NS equation in M'^ which going in Fourier space can be written 

7;t(fc) =e-l'=l'*^;o(fc)+^ Te-I^l'^*-"' / dk' {k,Vs{k ~ k'))PkVs{k') ds (33) 

where vt is the Fourier transform of the velocity field, (•, •) is the scalar product in C'^ and Pfe : C'^ — s- 
is the projection on the directions orthogonal to the vector k G M^, i.e. P^a = a — {k, a)fc|fc|~^. 
Eq. dSSl) will be studied in the spaces $(a), a e [2,3) where v € $(a) if w G C(R^;C^) with 
k-v{k) = and ||w||q, = supj.gjj3 |fc|"|?;(A;)| < oo. We will write a = 2 + e with e G [0, 1). The spaces 
$(a) can contain solutions with infinite energy and enstrophy so classical results about existence 
and uniqueness do not apply. Sinai [El IHl [ITj , studied eq. ((33|) in <i>(a) with a > 2, showing that 
there is existence of unique local solutions and that these solutions survive for arbitrary large time 
if the initial condition is small enough. Related works on NS are those of Le Jan and Sznitman [TB] , 
Cannone and Planchon |19j and the reviews of Bhattacharya et al. in and Waymire [H] . 

Following [6] we will describe the representation for these solutions as series indexed by planar 
binary trees. The use of trees to rigoroulsy analyze the NS equation has been somewhat pioneered 
by Gallavotti [22]. 
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The NS equation can be cast in the abstract form 

Ut = StUQ+ St-sB{us,Us)ds. (34) 
Jo 

where S" is a bounded semi-group on <i>(Q!) and _B is a symmetric bihnear operator which is usually 
defined only on a subspace of $(a). Here we cannot proceed as in the KdV case by going to the 
interaction picture since S is only a semi-group, so we must cope with the convolution directly. 
In [6] we showed that the solutions of this equation in the case of the 3d NS equation have the 
norm convergent series representation 

ut = S,uo+ ^toK'^^^) (35) 

tEBT 

where 0{t) is a degree function defined by 9{») = 2, 6l([r]) = 1 + 9{t), 0{[tiT2]) = e{T^) + e{T^) and 
the ^(r)-multilinear operators X'^ have recursive definition 

^tsif^"^) ^ / St-uB{Su-s^,Su-s^p)du 

J s 

xi:^^xw^^)+i)) = ^'5._„s(x::(^x«(^^)),^)d. 

and 

xr^^i(^><(«(^^H^(^^))) = ^*5,_„i?(x::(^><^(^^)),x::(^><^(-^)))du 

and by induction we can prove that, for any e G [0, 1) these operators are bounded by 

-|fc|2(t-s)/(|r| + l) 

\Xl{h''(^^^){k)\ < a {t - sp^/'WhfJ^^ (36) 

where the constants Cr can be chosen as 

= AI"I(t!)-"/2_ (37) 

for some other constant A > depending only on e. 

Due to the presence of the convolution integral these X operators does not behaves nicely with 
respect to the coboundary S. In 12 we introduced cochain complex (C*,^) adapted to the study 
of such convolution integrals where the coboundary S is obtained from (5 by a " twisting" involving 
the semigroup. There exists also a corresponding convolutional sewing map A which provide an 
appropriate inverse to 6. Algebraic relations for these iterated integrals have then by-now familiar 
expressions, e.g.: 

~SXl{^-^) = X'{X'{^-%Sy,) 

and so on. Indeed X can be considered as a branched rough path and the asymptotic behaviour (j37p 
supports somewhat the conjectured asymptotics (|20p . 

The series representation ([35]) together with the bounds (I36p imply that 

\ut{k)\ < \S,uo{k)\ + J2 ^T7^ i^'^'/'holl^^ (38) 
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By induction wc can prove that 7(t) > 21'^' ^ where equahty holds for the binary trees for which 
every path from the root to the leaves has the same length. This estimate together with 

(|r| + l)/2<0(r) < |r| + l (39) 

easily proven by induction on |t| give 

-(fe^t)/(n+l) 

\utik)\ < \StUo{k)\ + ^ Z„i3" — t-/2hollL"+^)/'(l + ||^^oIU)("+^)/' (40) 

n>l ' ' 

for some different constant B and where Z„ is the number of trees in BT with n vertices for which 
we have the estimate Z„ < D"{n + 1)~^/^ for a constant D > 0. So the series ([35]) is controlled by 
the geometric series PO]) and converges in norm t is small or \\h\\a is small. In the case e = the 
dependence in time of the r.h.s. is bounded and so the series converges for all time if the initial 
condition is small enough. The series gives also additional informations on the global solution when 
e = 0: 

a) dissipation: for fixed k G R'^\{0}, limi^oo |wt(fc)| = 0; 

b) smoothness: for fixed t > 0, there exists two constants C3, C4 such that \vt{k)\ < C'ae"'^*''^!^ 
as |fc| — > 00. 

In particular this second property can be proved by the Laplace method applied to the majorizing 
seriej40l 

Another interesting way to analyze the NS series ((35|) is to note that different classes of trees 
give different contributions. We define simple trees the trees with at most one branch at each 
vertex, i.e. of the form [•••[•]■••]• Short trees are instead trees for which at each vertex we have 
two branches, each of which carries (asymptotically) a fixed proportion (a or 1 — a) of the vertice 
and without loosing generality we consider a G (0, 1/2). We will denote BTq the set of simple trees 
and BTa the set of short trees corresponding to the proportion a. The distinction between these 
classes of trees is relevant when discussing the asymptotic behavior of the tree factorial. Indeed for 
r £ BTq we have t! = |t|! while for any a £ (0,1/2) there exists constants Di, D2, D3, D4 such 
that, for any r G BTa we have 

D3\t\-^d['U j{t) < D,\t\-^dI^K 

This different behavior is responsible for different convergence properties of the sum l|35p when 
restricted to simple or short trees. Indeed it is possible to prove that the series restricted to simple 
trees is convergent for all times whatever the size of the initial condition while the estimate for short 
trees do not ensure this important property. In some sense the difficulty of finding global solution 
of NS is due to the presence of arbitrarily large short trees in the expansion. A similar phenomenon 
is observed in [17J and exploited in [23] to prove blow-up for the complex Navier-Stokes equation 
by a renormalization group argument. 
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